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We investigate consequences of the effective colour-dielectric formulation of 
lattice gauge theory using the light-cone Hamiltonian formalism with a trans- 
verse lattice 0]. As a quantitative test of this approach, we have performed 
extensive analytic and numerical calculations for 2 -I- 1-dimensional pure gauge 
theory in the large N limit. We study the structure of coupling constant space 
for our effective potential by comparing with results available from conven- 
tional Euclidean lattice Monte Carlo simulations of this system. In particular, 
we calculate and measure the scaling behaviour of the entire low-lying glueball 
spectrum, glueball wavefunctions, string tension, asymptotic density of states, 
and deconfining temperature. 

The recent Euclidean Lattice Monte Carlo (ELMC) simulations of Teper |^ 
have shown that pure non-Abelian gauge theory behaves much the same way 
in three as in four dimensions: a discrete set of massive boundstates are gen- 
erated by a linearly confining string-like force. Moreover, Teper has performed 
calculations for = 2, 3, and 4, allowing an extrapolation to large N. The 
large-A^ limit is convenient, though not essential, for our Light-Cone Trans- 
verse Lattice (LCTL) formulation. Our formulation offers the rare possibility 
of describing the parton, constituent, and string behaviour of hadrons in one 
framework. The relationship between these pictures, each very different but 
equally successful, remains one of the outstanding enigmas of QCD. 

We characterise a dielectric formulation as one in which gluon fields, or rather 
the SU{N) group elements they generate, are replaced by collective variables 
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which represent an average over the fluctuations on short distance scales, rep- 
resented by complex N x N matrices M. These dielectric variables carry colour 
and form an effective gauge field theory with classical action minimised at zero 
field, meaning that colour flux is expelled from the vacuum at the classical 
level. The price one pays for starting with a simple vacuum structure is that 
the effective action will be largely unknown and must be investigated per se. 

Starting with the Wilson lattice action we take the continuum limit in 
the and directions, leaving the transverse direction discrete. Replacing 
the link variables U G SU{N) with N x N complex matrics M, one derives a 
transverse lattice action whose form was first suggested in Ref. Q 

I dx'^dx^ I] I^Tr {D^M,i{D"M,i)^} " ^ Tr {F.^pF'"'} - V^i [M] 

(1) 

where a,(3 £ {0, 2} and 

DaM^i = [da + iAa{x^)) M^i - iM^iAaix^ + a) . (2) 

M^i lies on the link between x^ and x^ +a while Aa(x^) is associated with the 
site x^. V^i [Af] is a purely transverse gauge invariant effective potential. Next 
we introduce light-cone co-ordinates x^ — x^ = [x'^ ± x^) / and quantise by 
treating x'^ as canonical time. The theory has a conserved current 



J^l — I 



i?" Ml + Aft D'^ A4i_, 



(3) 



If we pick the light-cone gauge A- = the 
non-propagating field A^ satisfies a simple constraint equation at each trans- 
verse site (d^)'^ A+{x^) — g'^J^i/a. Solving this constraint leaves an action in 
terms of the dynamical fields Af^i 

A^ [ dx+dx- J2 Tr |aaA4ia"A/t, + |! J+ j^J^i \ - V^i [M] . (4) 



At large N, Eguchi-Kawai reduction introduces considerable simplification. 
For = the theory is isomorphic to one compactified on a one-link transverse 
lattice, id est we can simply drop the argument x^ or I from M in all of the 
previous expressions. Effectively one is now dealing with a 1 -I- 1-dimensional 
gauge theory coupled to a complex scalar field in the adjoint representation 
(with self-interactions). 

For the transverse effective potential F[Af], we will include all Wilson loops 
and products of Wilson loops up to fourth order in link fields M: 

V[M] = ^2 rpj. l^^^/t} + ^ Tr {AfAftATAft} 

+ ^ Tr {AfMtA/tA/} + Tr {M^ M} Tr {Af^Af} (5) 
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Note that the last term above, which might appear suppressed at large N, is 
in fact non-zero only for 2 particle Fock states. 
Let us introduce creation/annihilation operators 




In the associated Fock space, we include only states annihilated by the charge 
/ dx~J^. This gives a Hilbert space formed from all possible closed Wilson 
loops of link modes a± on the transverse lattice. Thus, a typical p-link loop 
will be something like 

Tr{a^i(fci)aLi(fc2)ali(fc3)---aii(fcp)} |0) (7) 

where the number of +l's equals the number of — I's, J2m=i ^"i — P^, and 
km > 0. At large N we need only study the dynamics of single connected 
Wilson loops in the Hilbert space since the loop-loop coupling constant is of 
order 1/iV. These loops may be thought of as 'bare' glueballs, and the problem 
is to find the linear combinations that are on mass shell. Neglecting km = 0, 
which is consistent with expanding about the M = solution of the dielectric 
regime, the Fock vacuum is an eigenstate of the full light-cone Hamiltonian 
p- 10) = P+ |0> = 0. 

The theory possesses several discrete symmetries. Charge conjugation in- 
duces the symmetry C : a^^^ ^ ^-iji- There are two orthogonal reflection 
symmetries Vi and V2 either of which may be used as 'parity'. If Vi : x^ 
—x^, we have Vi : ij ^ '^-lij- rotational symmetry has been re- 
stored in the theory, states of spin 7^ should form degenerate Pi doublets 
I + ^7) ± I — J') We use "spectroscopic notation" to classify states. 

For Ai and A2 small there is very little mixing between Fock states of different 
number of link modes p. In this case a mass eigenstate |^) has predominantly 
a fixed p, the mass increasing with p. For a given p, the energy also tends to 
increase with the number of nodes in the wavefunction / due to the J (9-) ^ J 
term (0), which is in fact a positive contribution. Thus one expects the lowest 
two eigenstates to be approximately 

dfc/+i,_i(fc,P+-A:) Tr|a^i(fc)aLi(P+-fc)||0) (8) 

with the lowest state having a symmetric wavefunction — fc), cor- 

responding to 0'^'^, and first excited state having _i antisymmetric with 
one node, corresponding to 0~~. The next highest states are either a 4-link 
state with positive symmetric wavefunctions and or 

a symmetric 2-link state with having two nodes. In the glueball spec- 

trum we identify the latter states as 0;;!"+ and 2++, respectively, although actual 
eigenstates are a mixture of these. 
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Fig. 1. — Mass = jj,'^a/{g'^N) such that the lowest Af^ eigenvalue is zero vs 
\i/{g^N) and M/ig'^N) for p < 6, /S' = 14, A3 = lOOg^W. Below this surface the 
spectrum is tachyonic and below = our quantisation breaks down. The dark 
line is the scaling trajectory. 



In order to fix the coupling constants in the effective potential, we perform 
a least fit to Teper's ELMC large N extrapolated spectrum. As we shall 
later show, the mass in units of the coupling m? — jjp'a/ [g^N] is a measure of 
the lattice spacing while the other terms in our effective potential Ai, A2, and 
A3 must be found from the fitting procedure. We will also determine g^N/a 
based on a fit, which we check self-consistently with measurements of the string 
tension. 

In our numerical solutions we restrict the number of link fields p in our basis 
states and discretise momenta by demanding antiperiodicity of the fields in 
x~ — > + i. In order to minimise errors associated with these truncations, 
we take the spectra for various K and p truncations and extrapolate to the 
continuum, K,p ^ 00. Since we cannot measure \J'\ directly, we only classified 
states according to Vi and C during the fitting process. We estimate our one- 
sigma errors from finite K and p-truncation to be roughly 0.05M^. As we fit 
the various couplings as a function of , we find a narrow strip in parameter 
space where we obtain good agreement with the ELMC spectrum; see Fig. |l|. 
As we shall see, moving along this strip corresponds to changing the lattice 
spacing a. The strip, where has a local minimum, disappears when is 
sufficiently large, indicating that for large enough lattice spacing our truncation 
of the effective potential is no longer a good approximation. 

The numerical bound from absence of tachyons is shown in Fig. |l| as a zero- 
mass surface. As the transverse lattice spacing vanishes the mass gap should 
vanish in lattice units. The fixed point for this, which we believe lies somewhere 
at negative m?, should lie on the zero-mass surface, but is inaccessible to us 
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Fig. 2. — A fit of our LCTL extrapolated results against Teper's ELMC large A'' 
extrapolated spectrum. The eigenvalues are shown in units of (Teper's) string 



tension for various \ J\ 
„2 AT \ ^ _ 1 nn„2 



The parameters are m 

.2 / 



0.065, 

,2 



Ai 



0.2023'' A, Aa = 

0.55ff^Ar, As = lOOg^Af, and g^N = 3A7aa with = 23.5. One finds similar 
spectra along the entire scaling trajectory in coupling constant space. 



in the dielectric regime > 0. Nevertheless the scaling trajectory should 
gradually approach the zero-mass surface if it is to eventually encounter the 
fixed point. 

In Fig. 1^, we have plotted a typical spectrum along the scaling trajectory 
together with the ELMC results. For graphing purposes, we assigned \ J'\ to 
our spectrum based on a best fit to Teper's results. Although the overall fit 
with the conventional lattice results is quite good, we see two deficiencies of 
our spectrum that cannot be attributed to K or p-truncation errors. First, 
we see that the lowest state is too low in energy. Second, we see that 
the lowest parity doublet 2^~^ is not quite degenerate. We believe that these 
discrepancies must be due to our truncation of the effective potential. Finally, 
we have made no prediction for the lowest 1+'^ state since it lies too high in 
the IJ'I''"'' spectrum. 

To measure the string tension in the direction (before Eguchi-Kawai re- 
duction) consider a lattice with n transverse links and periodic boundary con- 
ditions. Constructing a basis of Polyakov loops, "winding modes," that wind 
once around this lattice, one may extract from the lowest eigenvalue vs 
n the lattice string tension aa — AM„/An. String theory arguments indi- 
cate that oscillations of the winding mode transverse to itself yield a form 
M^/cr = aa?"n? — 7r/3, the constant correction being due to the Casimir en- 
ergy |6|. Fig. ^ shows a typical vs n plot for winding modes where we see 
a good fit to a quadratic. The constant term -3.5, however, does not agree well 
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(a) (b) 

Fig. 3. — (a) Lowest eigenvalue vs n for winding modes. Here, K = 10.5 or 11, 
p < n + 4, and the couplings are taken from Fig. ^ Also shown is a numerical fit to 
1.92n^ — 3.5. The dashed line is from an analytic estimate, (b) The lattice spacing 
in units of the string tension along the scaling trajectory. Here we plot a^a vs m^. 



with the expected value of — 7r/3. 

As a consistency check, we use aa^ from the string tension measurement 
together with g'^N/ (acr) from the spectrum to form the ratio ag^ (as measured 
by us) to ag^ (as measured by ELMC). The result is equal to one ±5% all along 
the scaling trajectory. If we then assume that our a is equal to the ELMC value 
of (T, we can determine the lattice spacing a in units of the string tension. This 
demonstrates that the mass — ii^a/{g'^N) determines the lattice spacing 
and that the continuum limit occurs at < 0; see Fig. ||. 

Another quantity of interest is the deconfinement temperature Tc which may 
be obtained from the Hagedorn behaviour of the asymptotic density of 
mass eigenstates p{M) ~ M~°' exp{M/Tc). The canonical partition function 
diverges for T > Tc. li a > {D + l)/2 it is a phase transition, beyond which 
the canonical and microcanonical ensembles are inequivalent. If a < (D + l)/2 
the ensembles are equivalent and Tc represents a limiting temperature — the 
free energy diverges at Tc- It is essential to demonstrate that the spectrum 
is sufficiently converged in K, thus we only fit to the states between 0.5 < 
log(t) < 5.5 for the K — 7 data in Fig. ^. A numerical fit gives. 

Since the density of states is p{M) = dt/dM, we find Tc = Q.78y/a- with an 
estimated error of at least 10%. Due to the large error, this result is compatible 
with the Euclidean lattice result j|, ||. The fact that the power correction 
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a w 4.86 in the density p{M) is much larger than {D + l)/2 — 2 means 
we can safely say that the thermodynamic free energy remains finite at Tc, 
marking a true phase transition (presumably deconfinement) rather than a 
limiting temperature. In addition, we have an analytic estimate suggesting 
that Tc lies in the range 0.81y/a < Tc < 0.98-^^, also in agreement with the 
ELMC result. 



References 



[1] S. Dalley and B. van de Sand e, Cambridge University report No. 
DAMTP-97-16, bep-ph/9704408i 



[2] M. Teper, Phys. Lett. 289B (1992) 115; Phys. Lett. 311B (1993) 223; 
Nud. Phys. (Proc. Suppl.) B53 (1997) 715; Ox ford University Report 
OUTP-97-01P (unpubHshed) |hep-lat/9701004| . 
[3] K. G. Wilson, Phys. Rev. D 10 (1974) 2445. 

[4] W. A. Bardeen and R. B. Pearson, Phys. Rev. D 14 (1976) 547; W. A. 

Bardeen, R. B. Pearson, and E. Rabinovici, Phys. Rev. D 21 (1980) 1037. 
[5] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1983) 1063. 
[6] M. Liischer, Nucl. Phys. B180 (1981) 317; P. de Forcrand, H. Schneider, 

M. Teper, Phys. Lett. B160 (1985) 137. 
[7] R. Hagedorn, Nuovo Cimento Suppl. 3 (1965) 147; S. Frautchi, Phys. 

Rev. D 3 (1971) 2821; R. D. Carlitz, Phys. Rev. D 5 (1972) 3231. 
[8] M. Billo, M. Caselle, A. D' Adda, and S. Panzeri, Int. J. Mod. Phys. A 

12 (1997) 1783. 



